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ABSTRACT 

The flat spectrum radio quasar (FSRQ) PKS 2123—463 was associated in the First 
Fermi-LAT source catalog with the 7-ray source 1FGL J2126. 1—4603, but when con- 
sidering the full first two years of Fermi observations, no 7-ray source at a position 
consistent with this FSRQ was detected, and thus PKS 2123—463 was not reported in 
the Second Fermi-LAT source catalog. On 2011 December 14 a 7-ray source position- 
ally consistent with PKS 2123—463 was detected in flaring activity by Fermi-LAT . 
This activity triggered radio-to-X-ray observations by the Swift, GROND, ATCA, Ce- 
duna, and KAT-7 observatories. Results of the localization of the 7-ray source over 41 
months of Fermi-LAT operation are reported here in conjunction with the results of 
the analysis of radio, optical, UV and X-ray data collected soon after the 7-ray flare. 
The strict spatial association with the lower energy counterpart together with a simul- 
taneous increase of the activity in optical, UV, X-ray and 7-ray bands led to a firm 
identification of the 7-ray source with PKS 2123—463. A new photometric redshift 
has been estimated as z = 1.46 ± 0.05 using GROND and Swift/XJVOT observations, 
in rough agreement with the disputed spectroscopic redshift of z — 1.67. We fit the 
broadband spectral energy distribution with a synchrotron/external Compton model. 
We find that a thermal disk component is necessary to explain the optical/UV emis- 
sion detected by Swift/LSVOT. This disk has a luminosity of ~ 1.8 x 10 46 erg s _1 , and 
a fit to the disk emission assuming a Schwarzschild (i.e., nonrotating) black hole gives 
a mass of ~ 2 x 10 9 M Q . This is the first black hole mass estimate for this source. 

Key words: galaxies: active - galaxies: nuclei - quasars: general - quasars: individual 
(PKS 2123-463) - gamma rays 



1 INTRODUCTION 

Blazars constitute the most extreme subclass of Active 
Galactic Nuclei (AGN), characterized by the emission of 
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strong non-thermal radiation across the entire electromag- 
netic spectrum and in particular intense 7-ray emission 
above 100 MeV. The typical observational properties of 
blazars include irregular, rapid, high-amplitude variability, 
radio-core dominance, apparent super-luminal motion, a flat 
radio spectrum, and high and variable polarization at radio 
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and optical frequencies. These features are interpreted as 
resulting from the emission of high-energy particles acceler- 
ated within a relativistic jet closely aligned with our line of 
sight and launched in the vicinity of the supermassive black 
hole harboured by the activ e galaxy (|Blandford fc Reesl 
ll97Sl ; IUrrv fc Padovani Ill995l ). 

Since the advent of the Energetic Gamma Ray Exper- 
iment Telescope (EGRET) on the Compton Gamma-Ray 
Observatory (CGRO), blazars were known to dominate the 
extragalactic high-energy sky. However, EGRET did not 
pinpoint the location of many sources with sufficient pre- 
cision to enable astronomers to associate them with known 
objects, leaving the legacy of a large fraction of unidentified 
sources in 7 rays. The point-spread function and sensitivity 
of the Large Area Telescope (LAT) on board Fermi pro- 
vides an unprecedented angular resolution at high energies 
for localizing a large number of newly found 7-ray emitting 
sources. Correlated variability observed at different frequen- 
cies can give important information for the identification of 
a 7-ray source with its low-energy counterpart. 

On 2011 December 14 a 7-ray flare from a source 
positionally cons istent with PKS 2123—463 w as detected 
by Fermi-LAT (|Orienti fc D'Ammandol 1201 ll). triggerin; 
GROND and Swift follow -up observations ( Rau et al. 1201 
iD'Ammando et "all l201l] ) that confirmed contemporaneous 
activity in the optical/UV as well as marginally in X-rays. 

PKS 2123—463 is a bright radio quasar with a lumi- 
nosity at 1.4 GHz Li.4GHz= (1.5±0.2)xl0 28 W/Hz (assum- 
ing the redshift estimated in this paper, z — 1 .46±0.0fJ3: 
see Sect. [BJ. On the basis of its spectral index ct r ~0.4 
(S(z/) oc iy~ ar ) between 408 MHz and 4.8 GHz, it was in- 
cluded in the CRA TES catalog of flat spectrum objects 
ijHealev et al. 1 120071 ). ATCA observations performed almost 
simultaneously at 4.8, 8.6 and 20 GHz during the Australia 
Telescope 20 GHz (AT20G) survey indicated a flattening of 
the radio spectrum at high frequencies to a spectral index of 
a r ~ 0.2. Polarized emission has been detected only at 4.8 
GHz where t he polarization is abou t 4% of the total intensity 
flux density l|Massardi et al.l 12008). A Chandra observation 
of PKS 2123-463 in March 2004 has show n the presence 
of a jet-like extended structure in X-rays (Marshall ct al. 
l201ll ). The redshift of PKS 2123- 463 was reported to be 
z = 1.67 (|Savage fc Wrightlll98ll ) based on an objective- 
prism spectrum, subsequently questioned by I Jackson et all 
(|2002h because two possible redshifts (0.48 and 1.67) were 
given from observation of two lines in the spectrum, and 
the motivation for the exclusion of the smaller redshift was 
not provided. This obj ect was a member o f the pie-Fermi 
launch Roma-BZCAT (|Massaro et alj 12009) catalog listing 
candidate 7-ray blazar s detectable by Fer mi-LAT but not in 
the CGRaBS catalog l|Healev et al.H200Sl ). 

The flat spectrum radio quasar (FSRQ) PKS 
2123- 463 (R.A :21h2 6m30.7042s, Dec.:-46d05m47.892s, 
J2000; IE et al. I l2006l ) was associated in the First Fermi- 
LAT s ource catalog (1F GL, 2008 August 4 - 2009 
July 4; lAbdo et all l2010aT l with the 7-ray source 1FGL 
J2126. 1—4603, while no association with a 7-ray source was 



1 The corresponding luminosity distance for z = 1.46±0.05 is Al 
= 10.6±0.5 Gpc and 1 arcsec corresponds to a projected distance 
of 8.53±0.02 kpc. 



present in the Second Fermi- LAT source catal og (2FGL, 
2008 August 4 - 2010 August 4; lNolan et ai1l2012i ). although 
a 7-ray source, 2FGL J2125. 0-4632 with a radius of the 95% 
source location confidence region of 0.17°, at 0.52° from the 
radio position of PKS 2123—463 was reported. Taking into 
consideration the high variability of blazars, the flux of PKS 
2123—463 could have decreased in the second year of Fermi- 
LAT operation. In addition, in the construction of the 2FGL 
catalog, PKS 2123—463 was split int o more than one c an- 
didate source seed (see Section 4.2 of lNolan et al.ll2012l . for 
details). These two factors may have led to the lack of its 
association with a 7-ray source. 

In this paper we present the localization over 41 months 
of Fermi-LAT data of a 7-ray source associated with the 
FSRQ PKS 2123-463. The correlated variability observed 
in optical, UV, X-ray and 7 rays confirms the identification. 
In addition a new estimation of the redshift of the source by 
means of the fit of simultaneous GROND and Swift/WOT 
data collected soon after the 7-ray flare is presented. The pa- 
per is organized as follows: in Section 2 we report the LAT 
data analyses and results. In Sections 3 and 4 we report the 
results of the Swift and GROND data analysis, respectively. 
Radio data collected by the ATCA, Ceduna and KAT-7 tele- 
scopes are presented in Section 5. Section 6 presents an esti- 
mation of the redshift of the source. In Section 7 we discuss 
the modeling of the overall SED and draw our conclusions. 
Throughout the paper a A-CDM cosmology with Ho = 71 
km s _1 Mpc -1 , Qa = 0.73, and Q m = 0.27 is adopted. 



2 FERMI-LAT DATA: SELECTION AND 
ANALYSIS 

The Fermi-LAT is a 7-ray telescope operating from 20 MeV 
to > 300 GeV. It has a large peak effective area (~ 8000 cm 2 
for 1 GeV photons), an energy resolution of typically ~ 10%, 
and a field of view of about 2.4 sr with an angular resolution 
(68% containment angle) better than 1° for energies above 
1 GeV. Further deta ils about the Fermi-LAT are given by 
lAtwood et all i|2009f ). 

The Fermi-LAT data reported in this paper were col- 
lected over the first 41 months of Fermi operation, from 2008 
August 4 (MJD 54682) to 2012 January 4 (MJD 55930). 
During this time the LAT instrument operated almost en- 
tirely in survey mode. The spectral analysis was performed 
with the ScienceTools software package version v9r23pl. 
The Fermi-LAT data were extracted from a circular Region 
of Interest (Rol) with a 15° radius centred at the radio loca- 
tion of PKS 2123—463. Only events belonging to the 'Source' 
class were used. The time intervals when the rocking angle 
of the LAT was greater than 52° were rejected. In addition, 
a cut on the zenith angle (< 100°) was also applied to reduce 
contamination from the Earth limb 7 rays, which are pro- 
duced by cosmic rays interacting with the upper atmosphere. 
The spectral analyses (from which we derived spectral fits 
and photon fluxes) were performed with the post-launch in- 
strument response functions (IRFs) P7S0URCE_V6 using an 
unbinned maximum likelihood method implemented in the 
Science tool gtlike. 

The background model used to extract the 7-ray sig- 
nal includes a Galactic diffuse emission component and an 
isotropic component. The model that we adopted for the 
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Figure 1. 2-month integrated flux (E > 100 MeV) light curve of 
PKS 2123-463 obtained from 2008 August 4 to 2012 January 4 
(MJD 54682-55930). The last bin covers 1 month of observation. 
Arrows refer to 2-er upper limits on the source flux. Upper limits 
are computed when TS < 10. 
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Figure 2. Integrated flux (E > 100 MeV) light curve of PKS 
2123-463 obtained from 2011 December 4 to 2012 January 4 
(MJD 55899-55930) with 2-day or 4-day time bins. Arrows refer 
to 2-0" upper limits on the source flux. Upper limits are computed 
when TS < 10. 



Galactic component, the same as used for the 2FGL catalog, 
is given by the file gal_2yearp7v6_v0.fits, and the isotropic 
component, which is the sum of the extragalactic diffuse 
emission and the residual charged particle background, is 
parametrized by the file iso_p7v6source.tx10. The normal- 
izations of both components in the background model were 
allowed to vary freely during the spectral point fitting. 

We examine the significance of the 7-ray signal from 
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the sources using the Test Statistic (TS) based on the 
likelihood ratio test. The Test Statistic, denned as TS = 
2Alog(likelihood) between models with and without the 
source, is a measure of the probability of having a 7-ray 
source at the localization specified, which compares models 
whose parameters have been adjus ted to maximize th e like- 
lihood of the data given the model (|Mattox et al.fl l996). The 
source model used in gtlike includes all the point sources 
from the 2FGL that fall within 20° of PKS 2123-463. 
The spectra of those sources were parametrized by power- 
law functions, except for 2FGL J2056.2— 4715 for which we 
used a log-parabola, and 2FGL J2124.6-3357 and 2FGL 
J2241.7— 5236 for which we used an exponentially cut-off 
power-law in their spectral modeling as in the 2FGL cat- 
alog. We removed from the model the sources having TS 
< 25 and/or fluxes (0.1-100 GeV) below l.OxlO" 8 photons 
cm" 1 s" 1 over 41 months and repeated the fit. We tested 
whether two distinct 7-ray sources (one at the radio posi- 
tion of PKS 2123—463 and one at the 7-ray position of 2FGL 
J2125. 0—4632) are detected simultaneously by Fermi with 
TS ^ 25 over 41 months of observations. In this scenario, the 
fit yields a TS of 553 for PKS 2123-463 and a TS of 14 for 
2FGL J2125.0-4632. We thus conclude that PKS 2123-463 
and 2FGL J2125. 0—4632 are the same source and we main- 
tain only PKS 2123—463 in the model. Thus a final fitting 
procedure has been performed with all sources within 10° of 
PKS 2123—463 included with the normalization factors and 
the photon indices left as free parameters. For the sources 
located between 10° and 15° we kept the normalization and 
the photon index fixed to the values obtained in the previ- 
ous fitting procedure. The Rol model includes also sources 
falling between 15° and 20° from the target source, which 
can contribute to the total counts observed in the Rol due 
to the energy-dependent size of the point spread function of 
the instrument. For these additional sources, normalizations 
and indices were fixed to the values of the 2FGL catalog dur- 
ing all steps of the fitting procedure. 

The 7-ray point source localization using the gtf indsrc 
tool over the photons extracted during the period 2008 Au- 
gust 4 - 2012 January 4 results in R.A. = 321.609°, Dec. = 
-46.076° (J2000), with a 95% error circle radius of 0.047°, 
at an angular separation of 0.024° from the radio position of 
PKS 2123-463 (R.A. = 321.628°, Dec. = -46.097°, J2000). 
This implies a strict spatial association with PKS 2123—463. 
The fit with a power-law (PL) model, dN/dE cx (E/E )" r , 
to the data integrated over 41 months of Fermi operation 
in the 0.1-100 GeV energy range results in a TS = 557, 
with an integrated average flux of (4.76 ± 0.35) xlO -8 
photons cm -2 s _1 , and a photon index F = 2.45 ± 0.05 
(with Eo fixed to 670 MeV). In order to test for curva- 
ture in the 7-ray spectrum of PKS 2123—463 we used as 
an alternative spectral m odel the log-parabo l a (LP), dN/dE 
oc (E /E )- Q - /31og(E/Eo) (lLandau et all 1 19861 ; iMassaro et all 
[200J) , where the parameter a is the spectral slope at the en- 
ergy Eo and the parameter j3 measures the curvature around 
the peak. We fixed the reference energy Eo to 300 MeV. The 
fit with a log-parabola results in a TS = 554, with spectral 
parameters a = 2.35 ± 0.09, and /3 = 0.06 ± 0.04. Applying a 
likelihood ratio test to check the PL model (null hypothesis) 
against the LP model (alternative hypothesis), the power- 
law spectral model is favoured, indicating that no significant 
curvature was observed in the average 7-ray spectrum. 
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Figure \T\ shows the 7-ray light curve of the first 41 
months of Fermi observations of PKS 2123—463 built us- 
ing 2- month time bins, with the exception of the final (1- 
month) data point. For each time bin the photon index was 
frozen to the value resulting from the likelihood analysis 
over the entire period. If TS < 10 2-cr upper limits were 
calculated instead. The systematic uncertainty in the flux is 
energy dependent: it amounts to 10% at 100 MeV, decreas- 
ing to 5% at 560 MeV, a nd increasing to 10% above 10 GeV 
l|Ackermann et al J 12013 ). 

As shown in Fig. [TJPKS 2123—463 was in a low 7-ray 
state (0.1-100 GeV flux < 5xl0~ 8 ph cm" 2 s" 1 ) for the 
first 2.5 years. In particular the source was not detected 
with 2-month time bin during the period 2009 August- 
2010 May. A first increase of the 7-ray flux was observed 
in 2011 February-March, and subsequently a flaring activ- 
ity in 2011 December. We show a light curve focused on 
the period of high activity (2011 December 4-2012 Jan- 
uary 4; MJD 55899-55930) with 2-day or 4-day time bins 
(Fig. [2]). The peak of the emission was observed between 
December 13 15:43 UT and December 14 15:43 UT, with 
a flux of (128 ± 23) xl0~ 8 photons cm -2 s _1 in the 0.1- 
100 GeV energy range, a factor of ~ 25 higher with respect 
to the average 7-ray flux observed during 2008-2011. The 
corresponding observed isotropic 7-ray luminosity peak in 
the 0.1-100 GeV energy range is 8.9±0.8 x 10 48 erg s" 1 
(assuming z — 1.46±0.05, see Sect. O, comparable to the 
values r eached by the most powerful FSRQs in flarin g activ- 
ity (e.g. lOrienti et alj|20ia lAckermann et alj|201(f ). Leav- 
ing the photon index free to vary during the period of high 
activity and -Bo fixed to 670 MeV the fit results in a pho- 
ton index F = 2.26 ± 0.06, showing a moderate harder- 
when-brighter beha viour already observed in other FSRQs 
(|Abdo et al.ll2010bl ). Replacing in the same period the PL 
with a LP, fixing the reference energy Eo to 300 MeV, we 
obtain spectral parameters a = 2.09 ± 0.13, and f3 = 0.10 ± 
0.04. We used a likelihood ratio test to check the PL model 
(null hypothesis) against the LP model (alterna tive hypoth- 
esis). These values may be compared, following iNolan et al.l 
(2012), by evaluating the curvature Test Statistic TS clir ue 
= 2(log(like)ipTog(like)pL)=3 (~1.7-<r); thus no significant 
curvature was observed in the 7-ray spectrum also during 
the high activity period. 

Finally, using the gtsrcprob tool, which calculates the 
probabilities an event belongs to each of the sources in the 
model, we estimated that the highest-energy photon emitted 
by PKS 2123—463 was observed with a probability of 87% at 
distance of 0.35° from the source with an energy of 92 GeV. 
For a single high-energy photon (^ 1 GeV), sky location can 
only be determined within ~ 0.8° (95% confidence radius). 



3 SWIFT DATA: ANALYSIS AND RESULTS 



The Swift satellite l|Gehrels et al.l [20041 ) performed two ob- 
servations of PKS 2123—463 in 2011 December triggered 
by the 7-ray flare detected by Fermi-LAT. As a compari- 
son we also analysed two additional observations performed 
during 2010-2011. The observations were performed with 
all three on-board i nstruments: the X-ray Telescope (XRT; 
iBurrows et~ai1l2005l. 0.2-10.0 key), th e UltraViolet Optical 
Telescope fUVOT; lRoming et al.ll2005l . 170-600 nm) and the 



Burst Alert Telescope (BAT; lBarthelmv et al.ll2005l . 15-150 
keV). 

The hard X-ray flux of this source is below the sen- 
sitivity of the BAT instrument for the short exposure 
time of these observations. Moreover, the source was not 
resent in the Swift BA T 58-month hard X-ray catalog 
Baumgartner et alj|2010) and the 54-month Palermo BAT 



catalog IjCusumano et al 



2010). 



3.1 Swift/XRT 

The XRT data were processed with standard procedures 
(xrtpipeline vO.12.6), filtering,and screening criteria by 
using the HEAsof t package (v6.1ljj. The data were collected 
in photon counting mode in all observations, and only XRT 
event grades 0-12 were selected. The source count rate was 
low (< 0.5 cnt s _1 ); thus pile-up correction was not required. 
Source events were extracted from a circular region with a 
radius of 20 pixels (1 pixel ~ 2'.'36), while background events 
were extracted from a circular region with a radius of 50 pix- 
els away from the source region. Ancillary response files were 
generated with xrtmkarf , and account for different extrac- 
tion regions, vignetting and point spread function (PSF) cor- 
rections. We used the spectral redistribution matrices v013 
in the Calibration database maintained by HEASARC. 

Considering the low number of photons collected (< 
200 counts) the spectra were rebinned with a minimum 
of 1 count per bin and the Cash statistic |Cash I Il979l ) 
was applied. We fitted the spectrum with an absorbed 
power-law using th e photoelectric absorption model tbabs 
l| Wilms et alj|200ch . with a neutral h ydrogen column fixe d 
to its Galactic value (2.34 x 10 20 cm~ 2 : lKalberla et alj|200_] ). 
The fit results are reported in Table [T] and Fig. [4] 

A marginal increase (1.5 a) of the X-ray flux with re- 
spect to the previous Swift/XI{T observations was observed 
on 2011 December 15, soon after the 7-ray flare. The source 
remained at a similar flux level in X-rays on 2011 December 
19. 



3.2 Swift/UVOT 

UVOT UV/optical imaging was obtained during all four 
Swift observations of PKS 2123—463. The photometry was 
carried out on pipeline processed sky images downloaded 
from the S wift data centejfl following the standard UVOT 
procedure (|Poole et alj|2008r ). Source photometric measure- 
ments were extracted from the UVOT imaging data using 
the tool uvotmaghist (vl.l) with a circular source extrac- 
tion region that ranged from 3'.'5-5"radius to maximize the 
signal-to-noise. In order to remain compatible with the effec- 
tive area calibrations, whic h are based on 5" aperture pho- 
tometry (|Poole et alj|2008l ). an aperture correction was ap- 
plied where necessary. This correction was at maximum 5- 
6 % of the flux, depending on the filter. The UVOT pho- 
tometry is presented in Table and Fig. f3] Contempora- 
neous UVOT observations in 2011 December 15 found PKS 
2123—463 about 0.5 mag brighter in v-, u-, and uil-band, 0.7 
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Table 1. Log and fitting results of Swift/XRT observations of PKS 2123—463 using a power-law model with Nn fixed to Galactic 
absorption. 



MJD 


UT Date 


Exposure Time 


r 


Flux (0.3-10 keV) a 








[sec] 




[xlO -12 erg 


cm 2 s - 1 ] 




55467 


2010-09-28 


2712 


1.47±0.30 


1.31±0.46 




55745 


2011-07-03 


2912 


1.61±0.28 


1.20±0.32 




55910 


2011-12-15 


4947 


1.68±0.18 


1.90±0.34 




55914 


2011-12-19 


5521 


1.26±0.18 


1.83±0.27 








a : Observed flux. 








. UVOT Photometry of PKS 2123-463. 












UT Date 






AB Magnitude" 








uvw2 


uvm2 


uvwl 


u 


b 


V 


2010-09-28 18:29-23:30 


19.78 ±0.09 


19.34 ±0.10 


19.17 ±0.10 


18.66 ± 0.09 


18.61 ± 0.14 


18.59 ± 0.20 


2011-07-03 17:15-23:29 


19.83 ±0.20 


19.28 ±0.14 


18.75 ± 0.11 


18.47 ±0.12 


18.81 ± 0.25 


18.21 ± 0.26 


2011-12-15 20:45-01:25 


18.97 ±0.06 


18.58 ±0.06 


18.24 ±0.05 


17.95 ± 0.06 




17.64 ±0.06 


2011-12-19 04:20-04:39 


19.17 ±0.12 


18.80 ±0.12 


18.51 ± 0.09 


18.04 ±0.07 


17.77 ±0.07 


17.82 ± 0.15 


a : Corrected for Galactic foreground reddening of Eb-v = 


0.03 mae dSchleeel et al.lll998l 


). 



mag in m2-band, and about 0.8 mag in ui2-band compared 
to the UVOT observation performed on 2011 July 3. 

4 GROND DATA 

On 2011 December 18, at 01:21 UT, PKS 2123-463 was ob- 
served wit h the Gamma-ray O ptical/Near-Infrared Detector 
(GROND; iGreiner et al.ll2008l ) mounted on the MPG/ESO 
2.2-m telescope at La Silla, Chile. Preliminary results have 
been reported in iRau et ail i|201ll ). GROND is a 7-channel 
imager that observes in four optical and three near-IR 
channels simultaneously. The data were reduced and anal- 
ysed with the stand ard tools and methods described in 
iKriihler" et all |200Sl ). Calibration was performed against 
an SDSS stan dard star field (g'r'i' z') and against selected 
2MASS stars i|Skrutskie et al.ll2006l ) (JHK S ). This resulted 
in 1-cr accuracies of 0.05 mag (g'r'i'z', J, H), and 0.07 mag 
(Ks)- All magnitudes are correcte d for Galactic foreground 
extinction of E B -v = 0.030 mag (jSchlegel et al.lll99Sl ) and 
are summarised in Table [3] 



5 RADIO DATA 

PKS 2123-463 is monitored by the TAN AMI program 
(Tracking Active ga lactic Nuclei with Austral Milliarcsec- 
ond Interferometry; lOiha et al.l |2010| ) at a number of ra- 
dio frequencies and resolutions with the Australia Telescope 
Compact Array (ATCA) and the Ceduna facilities. ATCA 
is observing this source every few weeks with "snapshot" 
observations at frequencies from 5.5 GHz through 40 GHz 
where each frequency is the centre of a 2 GHz wide band 
and the fluxes are calibrated against the ATCA p rimary flux 
calibrator PKS 1934-638 (|Stevens et al.l 120121 '). Each flux 
density has a 1-a uncertainty of 0.01 Jy. The Ceduna radio 
telescope in South Australia is monitoring PKS 2123—463 at 
6.7 GHz. Each flux den sity has a 1-a uncertainty of 0.3 Jy 
l|McCulloch et al.ll2005l ). No sign of increased activity of the 



Table 4. ATCA (5.5, 9, 17, 19, 38, and 40 GHz) and Ceduna (6.7 
GHz) observations of PKS 2123-463 in 2011. 



UT Date 


Frequency (GHz) 


Flux (Jy) 


2011-05-17 


5.5 


0.86 




9.0 


0.77 




17.0 


0.67 




19.0 


0.64 




38.0 


0.45 




40.0 


0.44 


2011-10-14 


5.5 


0.86 




9.0 


0.76 




17.0 


0.60 




19.0 


0.56 




38.0 


0.42 




40.0 


0.41 


2011-11-08 


5.5 


0.85 




9.0 


0.71 




17.0 


0.57 




19.0 


0.53 




38.0 


0.43 




40.0 


0.43 


2011-12-19 


6.7 


0.81 


2012-02-18 


6.7 


0.78 



flux density was detected between 40 and 5.5 GHz before 
the 7-ray flare (see Fig. [3} . 

In addition observations of PKS 2123—463 were made 
after the 7-ray outburst using the KAT-7 array (the proto- 
type for MeerKAT) in the Karoo. Since the array is still un- 
dergoing commissioning tests some uncertainties remained 
about the absolute calibration scale. To minimize these un- 
certainties the observations were taken at similar times of 
days (hence Local Sidereal Time range) in late 2011 Decem- 
ber to early 2012 January, with 4-5 hr durations. The very 
first observation (December 21) failed because the online 
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Table 3. GROND Photometry of PKS 2123-463. 



UT Date 


9' 


AB Magnitude" 
T> %> z' 


J 


H 


K a 


2011-12-18 01:21 


17.78 ±0.05 


17.52 ±0.05 17.43 ±0.05 17.24 ± 0.05 


16.75 ± 0.06 


16.53 ± 0.06 


16.24 ±0.09 



a : Corrected for Galactic foreground reddening of Eg_v = 0.03 mag (Schlcgcl et al.|[l998h . 



fringe stopping was not working properly, but as the array 
is very small (longest baseline 200-m) and the frequency low 
enough (1.822 GHz) the delay tracking corrections could be 
done offline with only marginal loss in signal to noise ratio. 
Subsequent observations were done in this mode. There were 
also different problems with various receivers over this time 
so M^coverage and maps were not directly comparable. The 
bandpass and gain solutions (on PKS 1934—638) were good 
enough and checks against a nearby phase reference source 
(PKS 2134—470) gave sufficient confidence in the visibilities 
and phases to make simple maps of the inner 0.5° of the 
primary beam and verify that no sources stronger than 100 
mjy were present. 

Assuming a flux density for PKS 1934-638 of 13.6 Jy 
at 1.822 GHz (as interpolated from ATCA models) the mea- 
sured flux density of PKS 2123-463 is 1.13 Jy on 2011 De- 
cember 24, 26, 29, and 2012 January 7 and 13. The uncer- 
tainties on the absolute flux density scale is 5%. As a com- 
parison we extrapolate the flux density at 1.82 GHz using 
the values at 2.7, 0.8, and 0.4 GHz from the Parkes Cata- 
logue (1990) and the Molonglo Sky Survey (2008). Assuming 
the spectral index 0.4 derived between 2.7 and 0.4 GHz, we 
obtain a flux density at 1.82 GHz of 1.0 Jy, indicating that 
no significant variation was observed during the 7-ray flaring 
activity at this frequency. 




1.5 1 , , , , . . . 

0.2 0.5 1 2 

Wavelength (|im) 



Figure 3. UV-near-IR SED composed of Swift/WOT observa- 
tions obtained on 2011 December 19 and GROND data taken on 
2011 December 18. The simple power-law model (solid line) and 
the AGN/galaxy hybrid template (dashed line) suggest a photo- 
metric rcdshift of 2 ^ 1.45 (see text for details). 

models we find 2 = 1.46 ± 0.05 (x 2 = 23, P z = 99.2). This 
is in rough agree ment with the init i al spe ctroscopic redshift 
of z = 1.67 from lSavage fc WrightJ (|l98ll ). 



6 NEW PHOTOMETRIC REDSHIFT 
ESTIMATION 

The redshift of PKS 2123—463 has not been estab lished con- 
vincin gly yet. Following the method described in lRau et al.l 
(2012) and applied to a sample of 103 Fermi-LAT blazars, 
we combined the GROND photometry with the UVOT pho- 
tometry from 2011 December 19 in order to construct a 13- 
band spectral energy distribution (SED). In order to ac- 
count for source-intrinsic variability between the GROND 
and UVOT observing epochs, the GROND photometry was 
corrected by —0.05 mag in all bands (see iRau et al.l I2OI2I. 
for a detailed description of the method). The resulting 
SED, corrected for the Galactic foreg round reddening of 
Eb-v = 0.03 mag (Schlcgcl ct al. 1998), has been fit with 
a set of power-law models as well as with hybrid tem- 
plates buil t from normal galaxies and AGN ([Salvato et al.l 
2009l,l201lf) using the Le PHARE code (|Arnouts et al.ll 19991 ; 
Ilbert et al.ll2006t ). The best-fitting redshifts for both tem- 
plate libraries are in good agreement (see Fig. [3}. For the 
power-law model we obtain a 99% confidence redshift of 
z = 1. 37tg'.2i (X 2 = 17, Pz = 94.6 0) and for the hybrid 



7 SED MODELING AND CONCLUSIONS 

Beyond the excellent spatial association obtained with the 
41-month -Fermi-LAT data set, the most secure and dis- 
tinctive signature for firm identification of the 7-ray source 
detected by Fermi-LAT with the blazar PKS 2123-463 is 
the simultaneous increase observed in the 7-ray, X-ray, and 
optical-UV bands (see Fig. [4] Sections [2] and [3]). 

In order to investigate the physical properties of the 
source we have built a simultaneous SED of the flaring state 
of PKS 2123-463. The Fermi-LAT spectrum was built with 
data from observations centred on 2011 December 10 to 19 
(MJD 55905-55916). In addition we included in the SED 
the GROND and Swift (UVOT and XRT) data collected 
on 2011 December 18 and December 19, respectively. Here, 
the flux suppression, in particular in the UV bands, was 
corrected assuming a power-law spectral shape consistent 
with the optical-near-IR measurements (Fa oc A 0,8 ). The 
data from ATCA and KAT-7 collected on December 19 and 
24, respectively, provided information about the radio part 
of the spectrum. The flare centred on ~ MJD 55909 had 
a variability timescale of ~ 2 days, which constrains the 
size of the emitting region during the flare to K' b < 2.1 x 
10 16 (fc/20) cm. 



5 The best fit templates and redshifts were selected on the ba- 
sis of a simple \ 2 method. P z is the integral of the probability 
distribution function J f(z)dz at 2 p hot ± 0.1(1 + 2 p j lot ), which 



describes the probability that the redshift of a source is within 
0.1(1 + 2) of the best fit value. 
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Figure 4. Light curves of PKS 2123—463 collected in 7 rays by 
Fermi-LAT (0.1-100 GcV flux in units of 10~ 8 ph cm" 2 s" 1 ), in 
X-rays by Swift/XRT (flux in units of 10 12 erg cm -2 s _1 ), in UV 
(filters wl: empty squares, m2: empty triangles, w2: filled circles; 
flux in units of 10 12 erg cm -2 s" 1 ), and in optical (filters u: 
empty squares, b: empty triangles, v: filled circles; flux in units of 
10 13 erg cm -2 s" 1 ) by Swift/WOT, and radio by ATCA (empty 
squares: 40 GHz, empty triangles: 17 GHz, filled circle: 5.5 GHz) 
between 2010 June and 2011 December. 



A "blue bump" accretion disc component is clearly vis- 
ible in the optical/UV data. We modeled these data with 
a combination of a nonthermal sync hrotron component and 
a Sha kura-Sunyaev disc component (|Shakura and Sunvaev I 
1973). A fit to the disc component allows us to get a rough 
estimate of the black hole mass of Mbh « 2 x 1O 9 M0 . This is 
the first black hole mass estimate for this source. It is con- 
sistent with black hole estima tes for other high-z F SRQs, 
obtained in a similar way by iGhisellini et alj (1201 if). Our 
mass estimate for PKS 2123—463 follows Ghisellini et al.l 
and assumes that the innermost stable circular or- 



bit (ISCO) is Rdisc,isco = 6-Rg, as one would expect for a 
Schwarzschild (i.e., nonrotating) black hole. If the jet is pro - 
duced by the Blandf ord-Znaiek dBlandford fc ZnaiekHl977l 1 
or Blandford-Payne (|Blandford fc PavneHl982T T mechanisms 
this requires a nonzero black hole spin, and one expects 
Rdiscjsco to be smaller or larger, dep ending on whether th e 
spin is retrograde or prog rade (e.g.. iGarofalo et all I2OI0I) . 
Due to the uncertainty in spin, this mass estimate can be 
considered to have considerable uncertainty. 

We modelled the portion of the SED from X-rays to 7 
rays assuming emission from a relativistic jet with mecha- 
nisms of synchrotron self-Compton (SSC), and Compton- 
scattering of a dust torus external to the jet (EC-dust). 
The d escri ption of the m o del c an be found in iFinke et al.l 
( 2008) and iDermer et all (|2009l ). The synchrotron compo- 
nent considered is self-absorbed below ~ 10 11 Hz. Correla- 
tions of 7-ray and optical flares with radio light curves and 
rotations of optical polarisation angles in low-synchrotron- 
peaked blazars seem to indicate the 7-ray/optical emitting 
region is outside the broad line region (BLR), where the dust 



torus is the likely seed photon source (e.g.. iMarscher et al.l 
2010). From the disc luminosity obtained by the fit of the 
blue bump, Ldi SC = 1-8 x 10 46 erg s _1 , we estimate an asso- 
ciated BLR radius Rblr = 4.4 x 10 17 cm, based on the re- 
lation between disc luminosity and Rbl r determined from 
reverberation mappin g campaigns (e.g., iKaspi et al. 1 120051 ; 
IGhisellini et al.1 2008). To minimize the scattered BLR con- 
tribution, we placed the emitting region at r > Rblr- Here 
the primary seed photon source is the dust torus, which 
was simulated as a one-dimensional ring with radius Rdust 
aligned orthogonal to the jet, emitting as a blackbody with 
temperature Td us t and luminosity Ldust- 

The model fit to the broadband SED can be seen in 
Fig.[5]and the parameters can be found in Table [5] A descri 
tion of the parameters can be found in lDermer et al.l (200! 
The dust parameters were chosen so that R dust is roughly 
consi stent with the sublimation radius (|Nenkova et al.l 
2008). The EC-BLR component was calculated assuming the 
seed photons are from Ha and have a luminosity of O.lLdisc- 

The electron distribution used, a broken power-law with 
index pi — 2.0 below the break at ~f' brk and p2 = 3.8 for 
Ibrk < 7' j i s consistent with particles injected with index 
2.8, and emission taking p lace in the fast cooling regime (e.g., 
iBottcher fc Dermeijl2002l) . That is, particles are injected be- 
tween "f' brk and y'max, with a cooling electron Lorentz factor 



"icool 



3m e c 



AcartescUtot 



where u to t is the total energy density in the blob frame, 
which in the case of our model fit is dominated by the ex- 
ternal energy density. In this case 7 coo ; is associated with 
"fmin, since it is in the fast cooling regime. Also note that in 
this model fit the magnetic field and electrons are nearly in 
equipartition. Jet powers were calculated assuming a two- 
sided jet. 

The Compton dominance for PKS 2123—463, i.e., 
the ratio of the peak luminosities of the Compton and 
synchrotron components, is ~ 50, whi ch is a rather 
stand ard value for powerful blazars (e.g., Ghisellin i et al.l 
2011). A large disc luminosity was estimated from the 
UVOT data, as expected for a powerful FSRQ, with a 
Ldi SC /Lg t j t j=0.2 in agreement w ith the blazar divide pro- 
posed by IGhisellini et all (|2009T ) as a result of the changing 
of the accretion mode. 

Variability is common in 7-ray blazars and provides 
a powerful tool to associate them definitively with objects 
known at other wavelengths and to study the emission mech- 
anisms at work. The combination of deep and fairly uniform 
exposure over ~3 hr, very good angular resolution, and sta- 
ble response of the Fermi-LAT is producing the most sen- 
sitive, best-resolved survey of the 7-ray sky. On the other 
hand, those cases where there is a decrease in the activity, 
and thus of the significance of detection, can lead to a more 
complex identification process for a 7-ray source. When com- 
bined with simultaneous ground- and space-based multifre- 
quency observations, the Fermi-LAT achieves its full capa- 
bility for the identification of the 7-ray sources with coun- 
terparts at lower energies and the knowledge of their emis- 
sion processes, as reported here for the high-z, Compton- 
dominated FSRQ PKS 2123-463. 
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Table 5. Model parameters for the SED shown in Fig. [5] 



Redshift 




1.46 


i_J LlliS. 1_JUI l^HLZi JT diLlUl 




20 


Dopplcr Factor 


Sr. 


20 


Meignetic Field 


jj 


U.O V .1 


Variability" Timescfile 


f 


1.7x10 s s 


V U 1 1 1 U V lllti 1 LLj KJ 1 1J1U l_/ 




4 1x1 16 cm 


Jet Height 




l.OxlO 18 cm 


Low- Energy Electron Spectrcil Index 


Pi 


2.0 


High-Energy Electron Spectral Index 


P2 


■ >.o 


1 n i rvi 1 1 r~n TV.l r^cV vt~~ir\ T /~woTi"f ^ h r^T" 

IV11I11111U.111 ±_Jlt_-Hjl Ull LiUl L-IILZj 1 aL LUI 


' min 


4 Q 


Break Electron Lorentz Factor 




l.OxlO 3 


Maximum Electron Lorentz Factor 


^fmax 


1.0x10 s 


Disc Luminosity 


-^disc 


1.8X10 46 erg s" 1 


Black Hole Mass 


Mm 


2 x 10 9 M 


Accretion efficiency 


V 


1/12 


Gravitational Radius 


Rg 


1.76X10 14 cm 


Inner Disc Radius 


^disc_ISCO 


6R 9 


Outer Disc Radius 




10 4 R 9 


Dust Torus luminosity 


Ldust 


1.5X10 45 erg s" 1 


Dust Torus temperature 


^dust 


1.7xl0 3 K 


Dust Torus radius 


-Rdust 


3.2xl0 18 cm 


Jet Power in Magnetic Field 




3.3X10 45 erg s" 1 


Jet Power in Electrons 


p 

1 j ,par 


1.5X10 45 erg s" 1 
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Figure 5. Spectral energy distribution data (circles and squares) 
and model fit (solid curve) of PKS 2123—463 with the model 
components shown as dashed curves. The data points were col- 
lected by GROND (2011 December 18), Swift (UVOT and XRT; 
2011 December 19), and Fermi-LAT (2011 December 10-19), to- 
gether with radio data from ATCA (2011 December 19) and KAT- 
7 (2011 December 24). 
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